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Current-driven magnetization dynamics of single ferromagnets in heavy metal/ferromagnet bi- 
layers has been recently realized. In this work, spin torque induced by so-called Rashba spin-orbit 
coupling and spin Hall effect are considered within a diffusive model. The dependence of the re- 
sulting torque as a function of the thicknesses of the ferromagnet and heavy metal is analyzed. We 
show that (i) both torques are on the form T = Turn x y + Tim x (y x m), (ii) the ratio Tii/T± 
strongly depends on the thickness of the layers and (iii) the thickness dependence of the spin torque 
provides an indication of the origin of the (Rashba- or spin Hall effect-induced) spin torque . 



PACS numbers: 72.25.-b,75.70.Tj,75.60.Jk 

I. INTRODUCTION 

The electrical manipulation of magnetic ultrathin lay- 
ers is presently attracting increasing interest due to the 
possible technological applications^^. Current-driven 
magnetization switching and excitations have been re- 
alized in spin-valves nanopillars and non-local configura- 
tion as well as magnetic domain walliP^. Interestingly, 
recent progress has been achieved towards the manipula- 
tion of single homogeneous ferromagnets in the absence 
of a polarizer or magnetic textur e 6 " 15 !. These new con- 
figurations take advantage of spin-orbit coupling to ma- 
nipulate the magnetization. An interesting example is 
the voltage-controlled anisotropy^, where the perpendic- 
ular magnetic anisotropy is controlled by an external gate 
voltage applied across a dielectric layer. In this config- 
uration, no current flow is involved and the microscopic 
origin of the anisotropy change is attributed to voltage 
control of the band filling^, also modeled by an effective 
interfacial Rashba spin-orbit coupling^. 

Alternatively, current-driven magnetization switching 
and excitation in single ferromagnets have been realized 
in the current-in-plane (CIP) configuration. These struc- 
tures consist of HM/F/MOx stacks where HM designates 
a heavy metal (Pt, Ta), F is a transition metal ferromag- 
net and MOx is a metal oxide (MgO or AlOx) (see Refs. 
I7 HT2")) . This effect has also been reported in strained 
dilute semiconductors^^. The observation of current- 
driven magnetization control in the absence of a spin po- 
larizer has been interpreted along two different schemes. 
One relies on the Rashba torque previously proposed in 
Refs. [TH] and [TH] and recently extended in Refs. [20J and 
I2T1 This torque, derived in the case of ferromagnetic 
2-dimensional electron gas (2DEG) in the presence of 
Rashba spin-orbit coupling^ is dominated by a field-like 
term along m x y (y is the in-plane direction transverse 
to the current injection and m is the magnetization direc- 
tion) and possesses a correction on the form of an (anti- 
)damping term along m x (y x m). Although applying 
the Rashba model to realistic trilayer structures involving 
ultrathin metallic layers remains questionable, Rashba 
spin-orbit coupling has been observed in a wide variety 



of metallic interfaces^!. Therefore, structures with asym- 
metric interfaces such as HM/F/MOx are appropriate to 
display interfacial Rashba interactions. 

Another possible mechanism relics on the Spin Hall 
effect (SHE) torque proposed by Liu et al^. In this 
case, the strong SHE present in the HM bottom layer 
injects a spin current normally to the interfaces with a 
polarization along j x z = y, where j is the direction of 
the current injection and z is the normal to the inter- 
faces. This SHE torque is equivalent to a torque arising 
from a polarizing layer that would be located below the 
F layer with its magnetization along y and the current 
injection along z. By pursuing the analogy, this torque is 
dominated by a (anti-)damping torque along m x (y x m) 
and a field-like torque along m x y. This field-like torque 
is usually disreg arded in spin- valves due to the short spin 
dephasing lengtfpSESl. However, when the layer thickness 
is comparable to the dephasing length, the spin current 
is not totally absorbed in the magnetic free layer and one 
can reasonably expect a field-like component to emerge^. 

Consequently, although the two torques, SHE and 
Rashba-driven, are expected to be dominated by different 
components (anti-damping torque for SHE and field-like 
torque for Rashba), they both adopt the same geometri- 
cal form T = T]|m x y + Tj_m x (y x m). A major issue is 
then to find a way to distinguish between the two origins. 
In this work, we model the spin transport in a HM/F bi- 
layer stack in three cases. In the first case, only SHE 
is allowed in the HM layer. In the second case, Rashba 
spin-orbit coupling is assumed to be present at the top 
F/MOx interface whereas SHE is absent in the HM layer. 
In the third case, Rashba is present at the bottom HM/F 
interface and SHE is neglected in HM. We find that the 
spin dynamics is different is these three cases leading to 
different thickness dependencies of the spin torque. 



II. TRANSPORT FORMALISM 

The transport formalism itself is subjected to discus- 
sion. The configuration of the samples under inves- 
tigation are thin multilayers in which a current is in- 
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jected in the plane (CIP). CIP configurations have been 
intensively studied both theoretically and experimen- 
tally in the early 90's, soon after the discovery of Gi- 
ant MagnetoResistance^ (GMR). The tremendous suc- 
cess of spin-valves configuration in which the current is 
injected perpendicular to the planed (CPP) has par- 
tially occulted the important discoveries established in 
CIP samples. From the theory viewpoint, whereas the 
semi-classical diffusion formalism developed by Valet and 
FerlpD (VF) has been tremendously successful in describ- 
ing CPP spin-valves, it simply does not apply in the 
case of CIP multilayers. One striking difference is that 
while spin transport is governed by the spin diffusion 
length A s / in CPP spin-valves, it is characterized by the 
spin-dependent mean free path \ a e in CIP multilayers. 
An instructive consequence of this seminal difference is 
that no CIP-GMR effect can be captured when apply- 
ing VF approach to CIP systems. The reason is that 
Bolztmann formalism, such as the one proposed by Cam- 
ley and BarnaJ^ (CB), treats the evolution of the spin- 
dependent carrier distribution function g a (t, r, v), which 
depends on the spin a, time t, position r and particle 
velocity v. In contrast, VF theory tracks the behavior of 
semi-classical quantities such as spin accumulation and 
spin current which implies averaging the carrier distri- 
bution over the velocity v. Therefore, whereas CB for- 
malism explicitly accounts for the particle velocity v, it 
disappears in VF theory (for details, see Ref. l2"§|) . 

To demonstrate the added value of Boltzmann for- 
mulation of the spin transport, we model a normal 
metal/ferromagnet bilayer in CIP configuration follow- 
ing Camley and Barna*^. In this approach, the semi- 
classical non-equilibrium electron distribution g a (t,r,v) 
is governed by Boltzmann transport equation 



dt 



dg a 



V r g a +v-V v g° = -■- 



(1) 



where r CT is the spin-dependent momentum scattering 
time. In a bilayer homogeneous along the (x-y) plane 
and perpendicular to the z direction, the solution of Eq. 
Q is on the form g a ± = ^^(1 + Ae ±z / T ^ v ^). This ap- 
proach has the advantage to explicitly account for the 
momentum distribution. In order to illustrate its im- 
plication on CIP transport, we solve Eq. ([I]) following 
Ref. 3TJ1 assuming diffusive interfaces and totally re- 
flective outer boundaries. The results are displayed in 
Fig. [T] together with the expected values in the diffusive 
limit, a — ne 2 T/m (dotted lines). Interestingly, where 
the electron mean free path A e = vpT is much smaller 
than the layers thickness (X e /d <C 1), the conductivity 
is very similar to the one expected in the diffusive limit. 
One the other hand, when the mean free path becomes 
on the same order as or larger than the layers thickness 
(X e /d 3> 1), the conductivity obtained within CB ap- 
proach strongly differs from the diffusive limit. 

We retain two important aspects from the results in 
Fig. [TJ (i) the current distribution is not homogeneous 
inside the layers and (ii) the effective conductance of the 
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FIG. 1. (Color online) Distribution of the conductivity in the 
direction normal to the plane of the bilayer. The dotted lines 
are the expected values in the diffusive approximation. The 
parameters are the same as in Ref. I3UI A^= 5 nm and A^,= 
0.5 nm. 



ultrathin layers is likely to be very different from the bulk 
conductances. Nevertheless, in realistic ultrathin layers 
the grains size and interfaces roughness are expected to 
play an important role and to reduce the effective mean 
free path. Therefore, item (ii) is not critical since the 
actual conductivity can be captured through an effective 
parameter. 

More importantly, the smooth variation of the current 
distribution at the vicinity of the interface may have 
quantitative consequences on the actual spin dynamics 
at the interface (where spin dephasing takes place^) and 
then on the magnitude of injected spin current. However, 
it can be argued that the conductivity jump at the inter- 
face remains quite sharp and therefore, the added value 
of Boltzmann formalism to accurately treat spin transfer 
torque in such bilayers may actually be limited. 

As a consequence, in the remainder of the paper we use 
a semi-classical diffusive model in line with VF theory. In 
general the spin and charge currents are written^ 

J = -VV <E> S - a H V [fx x Vn] <g) n (2) 
j e = -VVn-f3VV{S ■ m) + ala H VS x Vn (3) 

where J is the spin current density, S is the itinerant 
spin density (spin accumulation), n is the charge accu- 
mulation, T> is the diffusion coefficient, ajj is the Hall 
angle, /3 is the bulk spin asymmetry parameter and dp 
is the volume of the unit cell. The unit vector /i is the 
direction of the spin projection. The spin accumulation 
dynamics is given by 

dS _ _ 1 

dt 



= -V-J- 



1 s 

S x m m x (S x m) (4) 



where tj is the spin precession time, is the spin de- 
phasing length and t s j is the spin diffusion length. The 
system Eqs. Q-Q is solved assuming the following 
boundary conditions, illustrated in Fig. [2j 
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FIG. 2. (Color online) Schematics of the bilayer structure and 
boundary conditions for the three cases considered. 



• SHE only: total current reflection (^7=0) at both 
outer boundaries and continuity of spin current and 
spin accumulation at the HM/F interface. 

• Rashba at F/MOx interface: total current reflec- 
tion (i/=0) at the outer HM boundary, spin accu- 
mulation determined by Rashba-induced spin den- 
sity So at the F/MOx boundary and continuity of 
spin current and spin accumulation at the HM/F 
interface. 

• Rashba at HM/F interface: total current reflection 
{J=Q) at both outer boundaries and spin accumu- 
lation determined by Rashba-induced spin density 
S at the HM/F interface. 

We do not provide the detail of the derivation here. The 
spin accumulation is directly extracted from Eqs. Q-Q 
and the spin transfer torque is defined as 



. . . dSl\-V-J- 



-m 



(5) 



namely the spatial change of spin current, compensated 
by the spin relaxation term. f2 is the volume of the mag- 
netic layer. In the remaining of the article, we consider 
the magnetization oriented along z, i.e. out-of-plane. 



III. RESULTS 

A. General remarks 

Before entering in the details of the numerics, let us 
first comment on the general form of the spin transfer 



torque. The original concept of spin transfer torque pro- 
posed by Slonczewski relied on the transfer of transverse 
spin current to the local magnetization 1 , T = — V • J. 
This definition describes the core of the phenomenon, but 
in realistic structures, spin diffusion may modify the form 
of the spin torque. By exploiting the definition given in 
Eq. ^ and Eq. (g]), one can expression the resulting 
local torque in terms of spin current divergency 



T = -V -J 

Tj 



T sf \ T sf 



1 1 

1 — 

T <t> 



m x V • J, 



(6) 



Equation (|6| indicates that the spin torque is in general 
not simply given by the divergency of the spin current 
but possesses an additional component mxV-J arising 
from the presence of spin diffusion. Therefore, there is 
always a correction to the main torque component when 
spin relaxation is present. 

In the following sections, we address the dependence 
of the torques as a function of the layer thicknesses. The 
magnitude of the SHE torque is proportional to the cur- 
rent density in the heavy normal metal, ]hm, whereas 
the magnitude of the Rashba torque is proportional to 
the interfacial current jV In order to keep the contribu- 
tion of the spin dynamics explicit, we assume that jhm 
and ji are kept constant when varying the layers thick- 
nesses. In a realistic experiment, it is much easier to 
apply the same total current density. Assuming a sim- 
ple circuit description of the bilayer, one can relate these 
current densities to the total injected current density j't 



Jhm 
ji = 



3T, 



Pitt 

Oiti + ohm<1hm + (? F d F 



JT, 



(7) 
(8) 



where a v and d v (y = i,HM,F) are the conductivity 
and (effective) thickness of the interface, normal metal 
and ferromagnet, respectively. This renormalization pro- 
duces an additional thickness dependence that must be 
taken into account to properly describe the thickness de- 
pendence of the torque. Furthermore, in the light of the 
discussion about CB model of CIP transport, we also 
need to keep in mind that the effective conductivity of 
an ultrathin layer depends on the parameter X e /d, which 
implies that the conductivities cr„ are also thickness de- 
pendent. To avoid unnecessary complexity in the analy- 
sis of the results, we focus on the influence of the layers 
thickness on the spin dynamics itself (jN,i is kept con- 
stant). 



B. Case 1: Spin Hall effect Torque 

In this case, a spin Hall current polarized along y is 
generated in the HM layer and impinges onto the F layer. 
Fig. [3] displays the spin accumulation profile S x and S y 
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FIG. 3. (Color online) Case 1: Transverse spin accumula- 
tion S x (solid lines) and S y (dashed lines) profile along z in 
the presence of SHE in the heavy metal only, for different 
spin dephasing lengths. The other parameters are \ff — 15 
nm, Aj=0.5 nm, Af^ M = I nm, d,F~ I nm, cIhm=3 nm and 
V f /Vhm = f • 



FIG. 4. (Color online) Case 1: In-plane (T|| - solid lines) and 
out-of-plane torques (Tj_ - dashed lines) as a function of the 
F (a) and HM (c) layers thickness for different spin dephasing 
lengths A^>; spin torque ratio T±/T\\ as a function of the F (b) 
and HM (d) layers thickness. The parameters are the same 
as in Fig. [3] 



are a function of the position in the bilayer for different 
values of the spin dephasing. Following the definition of 
the spin torque in Eq. ([5]), the component S x produces 
the (anti-)damping torque and the component along S y 
produces the field-like torque. As expected, the S y com- 
ponent vanishes for very short spin dephasing. In this 
case, the transverse spin current J = x y®y is 
fully absorbed by at the HM/F interface and the torque 
reduces to the pure (anti-)damping torque. 

The thickness dependence of the resulting absolute 
torque c?T is represented in Fig. [4] Note that the torque 
is multiplied by the distance in order to the remove the 
1/d-dependence that naturally appears after volume av- 
eraging. As mentioned in the introduction, in the case of 
spin Hall effect, the in-plane torque (oc m x (y x m)) dom- 
inates over the out-of-plane torque (oc m x y). Whereas 
the absolute torque decreases for very thin F layers 



> 



I), it remains essentially unaffected by the 



Co thickness for dp > A^. 

Conversely, since the spin current is generated by 
asymmetric spin scattering in the bulk of the HM layer, 
increasing the HM layer thickness increases the spin 
torque in the range X^ M > d HM > 0. When the thick- 
ness of the heavy metal exceeds the spin diffusion length 
{dHM > Afj M ), the amount of spin Hall current injected 
into the ferromagnet saturates due to spin-flip scatter- 
ing (the additional HM thickness is inefficient in creating 
more spin Hall current). The thickness dependence is 
on the form oc 1 — cosh -1 duM /^f M ■, as proposed by 
Liu et a/P^. A last remark is that the ratio T±/T\\ is 
also thickness dependent, as shown in Fig. 4|b) and (d). 
While the ratio saturates towards the bulk value given 
byEq. ([§), it increases when decreasing the thickness of 
the ferromagnetic layer and decreases when decreasing 



the thickness of the heavy metal. Note however that a 
significant enhancement of the perpendicular component 
is only reached for extremely thin layers. Therefore, one 
can reasonably expect that the SHE torque reduces es- 
sentially to an in-plane torque (oc m x (y x m)) with 
a measurable dependence as a function of the HM layer 
thickness. 



C. Case 2: Rashba on the F surface 

In this case, the interfacial Rashba spin-orbit coupling 
is expected to produce current-induced non-equilibrium 



spin density along S = S^x + S^y 20 . For simplicity, 
we consider that the Rashba torque produces only 5° 
non-zero component (5° >> S® ~0). As shown in Fig. 
[5j this interfacial spin accumulation decays away from 
the surface producing a complex dynamics that generates 
both S x and S y spin density components in the ferro- 
magnet. For a vanishing spin dephasing length, only S x 
component survives which produces a field-like Rashba 
torque^. 

The thickness dependence of the resulting absolute 
torque dT is represented in Fig. [6| Whereas no sig- 
nificant thickness dependence is observed when varying 
the heavy metal thickness (Fig. |6jc)), a sizable thickness 
dependence is observed as a function of the thickness of 
the ferromagnet (Fig. |6^a)). Interestingly, we note that 
the in-plane torque can be significantly larger than the 
bulk value in the case of ultrathin ferromagnetic layers, as 
shown by the ratio T||/Tl in Fig. |6|b). This additional 
in-plane component has been theoretically shown to have 
a significant impact on the current-driven domain wall 
motion in the presence of Rashba torque^. Therefore, 
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FIG. 5. (Color online) Case 2: Transverse spin accumulation 
S x (solid lines) and S y (dashed lines) profile along z in the 
presence of Rashba spin-orbit coupling at the F/MOx inter- 
face only, for different spin dephasing lengths. The parame- 
ters are the same as in Fig. [3] 



FIG. 7. (Color online) Case 3: Transverse spin accumulation 
Sx (solid lines) and S y (dashed lines) profile along z in the 
presence of Rashba spin-orbit coupling at the HM/F interface 
only, for different spin dephasing lengths. The parameters are 
the same as in Fig. [3] 
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FIG. 6. (Color online) Case 2: In-plane (Tji - solid lines) and 
out-of-plane torques (T± - dashed lines) as a function of the 
F (a) and HM (c) layers thickness for different spin dephasing 
lengths A^; spin torque ratio Tu/T± as a function of the F (b) 
and HM (d) layers thickness. The parameters are the same 
as in Fig. [3] 



FIG. 8. (Color online) Case 3: In-plane (Th - solid lines) and 
out-of-plane torques (T± - dashed lines) as a function of the 
F (a) and HM (c) layers thickness for different spin dephasing 
lengths A^; spin torque ratio T\\/T±_ as a function of the F (b) 
and HM (d) layers thickness. The parameters are the same 
as in Fig. [3] 



while no thickness dependence is expected from the HM 
layer, tuning the F layer thickness may strongly enhance 
the Rashba-induced in-plane torque and have dramatic 
effects on current-driven magnetization dynamics. 



D. Case 3: Rashba HM/F interface 

The Rashba spin-orbit coupling at HM/F interface 
produces a non-equilibrium spin density along So = 
S°x + S°y « S°j0 that decays in both F and HM lay- 
ers. Since the source of the spin accumulation is confined 



at the HM/F interface, the decay in HM layer is inde- 
pendent on the spin dephasing length (see Fig. |7b. The 
thickness dependence of the resulting absolute torque dT 
is represented in Fig. [5J Again, the HM layer has no 
influence on the torque magnitude and reducing the F 
layer thickness can lead to an enhancement of the in- 
plane torque compared to the perpendicular torque, as 
shown in Fig. |8jb). Note however that this variation is 
only present in layers with thicknesses smaller than the 
spin dephasing length (d < A^). Therefore, in realistic 
systems, no significant difference from the bulk value is 
expected. 
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E. Discussion 

From the calculations above, we can draw three main 
conclusions. First, both Rashba- and SHE-induced 
torques are on the form T = 2]imXy+Tj_mX (y xm), the 
relative magnitude of Tji and T± is strongly dependent 
on the microscopic mechanism such as spin precession, 
dephasing and relaxation in the bilayer. Second, due to 
this complex spin dynamics, the ratio T]|/T_l strongly 
depends on the thickness of the layers. Third, while 
the three cases are affected by the ferromagnetic layer 
thickness, only the SHE-induced torque is affected by 
the thickness of the heavy metal. This influence is con- 
strained to thicknesses smaller than the spin diffusion 
length of the HM layer since the torque magnitude 
saturates beyond this length. 

Note that other contributions of the thickness depen- 
dence have been disregarded at this point and have to be 
accounted to accurately reproduce experimental values. 
First, as mentioned above, the torque has been evalu- 
ated at constant current density in the HM layer {jhm 
in case 1) and at the interface (ji in cases 2 and 3). This 
current density should be replaced by its expression in 
Eq. ([7]) or Eq. (J8|. In addition, at such ultrathin thick- 
nesses the effective conductivity of the layer depends on 
the ratio A e /d, as shown in Fig. [I] Finally, changing the 
thickness of the layers probably modifies the properties 
of the HM/F interface, such as interfacial resistivity and 



magnetic anisotropy, complexifying the analysis. 

However, although analyzing the effect of thicknesses 
variation on the spin torque requires a good understand- 
ing of the materials growth and its magnetic implications, 
the present study indicates that varying the HM layer 
thickness over a small range (on the order of Afj M ) is 
sufficient to identify the physical origin of the spin-orbit- 
induced torque. 

IV. CONCLUSION 

Using a semi-classical drift-diffusion description of the 
spin transport in a bilayer, both Rashba torque and 
SHE torque have been considered. The spin dynamics 
in such ultrathin layers has been investigated and the 
nature of the spin torque in such systems has been iden- 
tified. We showed that (i) both torques are on the form 
T = Xj|m x y + Tim x (y x m), (ii) the ratio T^/T^ 
strongly depends on the thickness of the layers and (iii) 
the thickness dependence of the spin torque may pro- 
vide an indication of the origin of the torque (Rashba- or 
SHE-induced). 
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